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Superlight small bipolarons from realistic long-range Coulomb and Frohlich 

interactions 
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We report analytical and numerical results on the two-particle states of the polaronic t- J p model 
derived recently with realistic Coulomb and electron-phonon (Frohlich) interactions in doped polar 
insulators. Eigenstates and eigenvalues are calculated for two different geometries. Our results show 
that the ground state is a bipolaronic singlet, made up of two polarons. The bipolaron size increases 
with increasing ratio of the polaron hopping integral t to the exchange interaction J p but remains 
small compared to the system size in the whole range < t/J p < 1. Furthermore, the model exhibits 
a phase transition to a superconducting state with a critical temperature well in excess of 100K. In 
the range t/J p < 1, there are distinct charge and spin gaps opening in the density of states, specific 
heat, and magnetic susceptibility well above T c . 

PACS numbers: 71.38.-k 



With few exceptions^, it is widely believed that the 
conventional Bardecn-Cooper-Schrieffer (BCS) theory 
and its intermediate-coupling Eliashbcrg extension 2 - do 
not suffice to explain high temperature superconductiv- 
ity. On the contrary, there is growing understanding 
that the true origin of high-temperature superconduc- 
tivity should be found in a proper combination of the 
Coulomb repulsion with a significant electron-phonon in- 
teraction (EPI)2. The many-body theory of strongly cor- 
related electrons and phonons was originally developed 
with the on-site Hubbard repulsion and the short-range 
Holstcin EPI using analytical strong-coupling expansion^ 
and powerful numerical techniques^ in the framework 
of the Hubbard-Holstein and Holstein-t J models^. Also 
the many-body Coulomb-Frohlich model, which takes 
into account a finite range of realistic interactions, was 
proposed^ and studied analytically^ and numerically^ 
showing a rich phase diagram with a polaronic Fcrmi- 
liquid, superconductivity induced by mobile bipolarons 
and a charge-segregated phase. In these and many other 
studies^ both interactions were introduced as input pa- 
rameters not directly related to the material. 

Recently it has been shown that, in highly polariz- 
ablc ionic lattices, the bare long-range Coulomb and 
electron-phonon interactions almost negate each other 
giving rise to a novel physics described by the polaronic 
t-Jp modelii with a short-range polaronic spin-exchange 
J p of phononic origin, 
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H = -53*i j <5 . ./ctc j +2 Jp(m-n) ( S m 

(i) 

Here the sum over n 7^ m counts each pair once only, 
S m = (1/2) J2a.cr> 4na^<?<?' c ma> is the spin 1/2 operator 
(r are the Pauli matrices), i = (m, a) and j = (n, a') 
include both site (m, n) and spin (a, a') indices; tij is the 
polaron hopping integral while J p (m — n) > t represents 



the exchange interaction between polarons on different 
sites from a residual polaron-multiphonon interaction. It 
has been proposed that the t- J p Hamiltonian, Eq.([T|), has 
a high-T c superconducting ground state protected from 
clustering^. 

In this work we present numerical and analytical re- 
sults on the two-particle eigenstates of the polaronic t- J p 
model as the building blocks for high-temperature su- 
perconductivity. It is worth noting that there is a wide 
difference between ([1]) and the familiar t-J modeli 2 - de- 
rived from the repulsive Hubbard U Hamiltonian in the 
limit U 3> t omitting the so-called three-site hoppings 
and EPI. The latter model acts in a projected Hilbert 
space constrained to no double occupancy. On the con- 
trary t-Jp Hamiltonian, Eq.([T]) has no constraint on the 
on-site occupancy since the on-site Coulomb repulsion is 
negated by the Frohlich EPI. The hopping integral Uj 
leads to the coherent (bi)polaron band while the antifcr- 
romagnctic exchange J p bounds polarons into superlight 
inter-site bipolarons. Moreover, the sign "+" instead of 
"— " in the last density-density interaction term in ([1]) 
provides an effective repulsion between pairs preventing 
their clusteringi 3 -, while the repulsive t-J model favors a 
phase separation. 

Also different from any model proposed so far, all 
quantities in the polaronic t-J p Hamiltonian (JT|) arc 
defined through the material parameters, in particular 
tij = T(m — n) cxp[— <? 2 (m — n)] with 
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(2) 



and 



J p (m) = T\m)/2g\m)hw , (3) 

where k = eoo£o/(£o~ e oo) and V is the normalization vol- 
ume. Here the high-frequency, and the static, eo di- 
electric constants as well as the optical phonon frequency, 
wq and the bare hopping integrals in a rigid lattice, T(m) 
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Figure 1: (Color online) Probability to find two polarons on 
the nearest-neighbor sites: Pb p (squares), on more distant 
sites P up (circles) and on the same site D (triangles) in the 
ground state of the t — J p Hamiltonian for chain (a) and zig- 
zag ladder (b). 
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Figure 2: (Color online) Two-particle ground-state energy Eo 
as a function of the hopping. Symbols correspond to ED 
(squares and triangles) and variational (crosses) data on finite 
clusters. We also report (circles) the results at fc = obtained 
by diagonalizing H(k) given in Eq.([S|. 



lowing two-particle singlet basis for an infinite lattice: 



are measured and/or found using first-principle Density 
Functional Theory-i^ in a parent polar insulator. 

In the following we restrict the range of the exchange 
interaction and of the hopping to nearest neighbors. One 
can readily fincU^ highly degenerate two-particle energy 
levels of the model ^ for t = 0, 

E (t = 0) = -J p , Ei{t = 0) = , E 2 (t = 0) = J p . (4) 

The ground and the highest energy states are bipolaronic 
spin-singlet and spin-triplet, respectively, made up of two 
polarons on neighboring sites. The zero-energy states are 
combinations of pairs of polarons separated by more than 
one lattice parameter and on-site bipolarons, since there 
are no on-site interaction terms in the Hamiltonian ([1]). 

For t ^ there is a finite bandwidth associated with 
each of the three energy levels. Exact diagonalization 
(ED) results show that the ground state configuration 
is virtually unchanged since the distance between two 
bound polarons remains of the order of the lattice spacing 
as long as t < J p . In that range in fact, regardless of 
the particular geometry, the probability Pb p to find two 
polarons on nearest neighbor sites decreases gradually 
with increasing t/J p and remains finite as shown in Figjl] 

We also show in FigJ2] the ground state energy as a 
function of the t/J p ratio for each analyzed geometry. 
Importantly, fitting the ED results in the t <C J P range, 
there is a contribution linear in t in the zig-zag ladder 
where a single hopping is sufficient for the coherent prop- 
agation of the intersite bipolaron through the lattice^. 
On the contrary, in the case of a one-dimensional chain 
the bipolaron hopping is realized through a second order 
process resulting in the quadratic behavior of the ground 
state energy as in the case of the on-site bipolaron^. As 
shown in Fig(21 ED results are in excellent agreement 
with the ones that can be obtained by means of the vari- 
ational method developed by Bonca et alJ^. 

Beyond ED results, additional information on the two- 
particle dynamic can be obtained by considering the fol- 



|m,k) = 



7n ?' 



(5) 

For the analyzed one dimensional geometries the matrix 
representation of the t-J p Hamiltonian in this basis is: 



( V2ei(fe) V2e 2 (fc) . . . \ 

\/2ei(k) -Jj, + V2e 2 (fe) ei(fe) e 2 (fc) '' 

V2e 2 {k) ei(fe) e 3 ei(fc) 

e 2 (fc) ei(fc) 



H(k) 



V 



(6) 



) 



where e\(k) = — It cos (fca/2), e<x{l£) = e% = for the 
chain, while we have ei(fe) = — 2t cos (fca/4), e2(fc) = 
—2t cos (fca/2), e3 = — J p for the zig-zag ladder. 

The eigenvalues of the tridiagonal matrix H (fc) deter- 
mine the energy dispersion E(k). In the limit t — > 0, 
ei(k) — > 0, so that the ground state energy is — J p in 
agreement with the ED results showed so far. For any 
t > 0, the problem is still solvable by requiring the sys- 
tem wave function to decay exponentially in the region 
where the potential vanishes^. In the case of a chain, the 
energy dispersion can be derived from a cubic equation 
that, in the t J p limit, gives: 



E a (k) = -Jp - (12t 2 /J p ) cos 2 (fca/2) + O (t 4 ) 



(7) 



with a quadratic contribution with respect to the hopping 
term. On the contrary, in the same limit the correspond- 
ing dispersion for the zig-zag ladder— has been found to 
be linear in t: 

E s (k) = - J p - t [cos (fca/2) + i/l + 4 cos 4 (fca/4)] + O (t 2 ) 

(8) 

As shown in Figj2j the energy dispersions at fc = ob- 
tained for the chain and the zig-zag ladder are in perfect 



3 



agreement with ED and variational results on finite clus- 
ters in the whole range <t/J p < 1. 

These results allow for some insight into a possible su- 
perconducting phase transition and pscudogap signatures 
in the response functions of the model. 

According to the Mermin Wagner theorem^, there 
should be no phase transition at finite temperatures in 
ID and 2D since there is no continuous symmetry break- 
ing. However, a finite temperature phase transition in 2D 
can exist via the Bcrezinsky-Kosterlitz-Thoulcss (BKT) 
mechanis m 18,19 . For a hard-core 2D Bose gas, where the 
Bose-Einstein condensation does not occu r 20 i 21 , a phase 
transition to a superfluid state is expected^. In particu- 
lar, it has been shownS 2 - in the dilute limit In ln(l/nf,r 2 ) > 
1 that the critical temperature is 

T c = 2irh 2 n b /k B m" lnln(l /n b r 2 ), (9) 

where rib is the boson density per unit area and r is the 
range of the boson-boson repulsion. In our case r ~ a 
and m** is the bipolaron mass. 

To estimate the bipolaron effective mass for a ID 
chain one can use our dispersion ((?]) that gives m** — 
h 2 J p /6a 2 t 2 . In the case of a zig-zag ladder, taking 
into account only the linear contribution in t/J p re- 
sults in an overestimated m** with m** /m* w 5 where 
m* = 2ft 2 /5ta 2 is the polaron mass^. Numerical results 
obtained by solving the eigenvalue problem of (|6]) show in 
fact that the ratio m**/m* is below this estimate, FigJ3] 
In particular, for t « J p we have m** ~ 2m* for both 
chain and zig-zag ladder geometries. 

Since the effective mass is proportional to l/t 2 or 1/t, 
one may conclude that T c should increase as t or t 2 with 
the polaron hopping integral. On the other hand, our ED 
results, Fig[TJ show that the probability Pb P {t) to find a 
hard-core tightly bound singlet decreases as the hopping 
increases. At low enough density BEC should not de- 
pend on whether the bipolarons arc nearest-neighbor or 
next-nearest neighbor, so long as they are bound and the 
bipolaron spacing is much greater than the typical po- 
laron separation. On the other hand when bipolarons 
overlap, their condensation appears in the form of the 
Cooper pairs in the momentum space with a lower crit- 
ical temperature, rather than in real space (BEC-BCS 
crossover—). Hence, bounds for the critical temperature 
can be estimated by weig hting Eq.© with P = P bp (t/J p ) 
as Tl rj Pb P {t/J p )T c . As shown in FigGH despite the low 
carrier density, the critical temperature is about 200K 
with the chain and the zig-zag ladder effective mass, for 
lib = 0.01/a 2 , a = 0.4nm and ftwo = 80meV. In particu- 
lar, obtained for the chain should not be considered 
as strictly related to the geometry but to the values of 
m** that could be a crude (but quite reliable) estima- 
tion of the bipolaron effective mass for a 2D lattice in 
the low-density limit. In the case of cuprate supercon- 
ductors with the polaron binding energy in the range 
0.5eV< E p < l.OeV and 0.3eV< J p < l.OeV -ii one 
gets the realistic value of the bare hopping integral 0.2 
cV< T(a) < 0.4cV that gives 0.05 < t/J p < 0.27 and the 
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Figure 3: (Color online) Ratio of bipolaron to polaron mass 
(left panel) in the t- J p model for different lattices and the re- 
sulting critical temperature (right panel) estimated with this 
mass at n;, = 0.01/a 2 and J p = l.OeV. For the zig-zag lad- 
der we report both the results obtained by using the linear-t 
dispersion— (dashed line) and the complete one (squares) cal- 
culated numerically from © by diagonalizing H(k). 



critical temperature 20A' <T C < 100/\ at rib = 0.01. 

Let us finally analyse the (pseudo)gap features in the 
density of states (DOS) and the spin susceptibility Xs 
of the polaronic t-J p model at high temperatures well 
above T£, when all carriers are non-degenerate. It is con- 
venient to introduce the "occupation density of states" 
(ODOS), p(uj, T) by weighting the standard temperature- 
independent DOS with the Fermi-Dirac and the Bose- 
Einstein distribution functions, 

p(w,T) = f.(u,T)M a (u) + 2f p (u,T)rf p (w) , (10) 

where: 

M» = ^-f ' dkS (lj - E 3 . P (k)) , (11) 

2t J-n/a 

and f s , p (uj,T) = [exp ((cj - fi s , p )/k B T) =F Here 
E SyP {k) is the (bi)polaron dispersion, fj, s = 2fx, fi p = [i are 
the chemical potentials of bipolarons and single polarons, 
respectively, with fj, < E s (0)/2. 

In the t = limit, according to ((H) we have three dif- 
ferent two-particle energy levels with E s (k) and E p {k) 
separated by J p . However, at low carrier density and 
temperature the highest energy level does not contribute 
to ODOS and we observe two sharp peaks at uj/J p = — 1 
and ui/J p = 0.0 with a suppression of ODOS around 
uj/J p — —0.5. Hence there is a single charge/spin pseu- 
dogap, A c = A s = J p . The ODOS for finite values of 
t/J p and temperatures is shown in Fig 2] with some gaus- 
sian broadening in the 8- function in Ea.(|ll[). modeling 
for instance disorder effects. At any t ^ 0, the two peaks 
become wider as t/J p increases and the gap between 
the bipolaron and the unpaired polaron bands gradually 
closes. With increasing temperature the single particle 
polaron band is more and more populated along with the 
incresing population of higher energy levels in the bipo- 
laron band, so that ODOS reflects a competition between 
bound and unbound states in the response functions. In 
particular, the behavior of the spin susceptibility indi- 
cates the presence of a finite spin-gap A s that decreases 
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Figure 4: (Color online) Signatures of a pseudogap opening in 
ODOS for different values of the polaron hopping (left panel) 
and temperature (right panel) calculated for the chain with 
a gaussian broadening 5 = 0.01J P , modeling a disorder effect 
in the 5-function in Eq. (|f f p . 



gradually as t becomes comparable with the binding en- 
ergy Jjr^, coherently with the suppression of the proba- 
bility to find a bipolaron in the nearest-neighbor configu- 
ration (FigJT]) . With increasing temperature, the number 
of occupied states within the pseudogap also increases. 
Without many-body correlation effects (i.e. a screening 
of the Coulomb and electron-phonon interactions at finite 
carrier densities^), the pseudogap itself does not depend 
on temperature; it is a matter of whether the temperature 
is high enough for single polaron states to have significant 
occupation. Hence a characteristic pseudogap tempera- 
ture T* exists in our model above which the pseudogap 



is suppressed, but it is a crossover temperature rather 
than a critical temperature. Further signatures of pseu- 
dogap opening are also found in the specific heatp 1 ^ where 
the Schottky anomaly is induced by thermal excitation 
within the bipolaronic and the polaronic bands. 

In conclusion, we have described some key features of 
the t — J p Hamiltonian in the low density limit. We have 
shown that the ground state configuration is a small bipo- 
laron singlet. Depending on the competition between the 
hopping t and the polaronic exchange interaction J p , the 
bipolaron size changes but remains small compared to 
the system size in the whole range < t/J p < 1. We 
have also argued that, in the 2D case, the presence of 
small light bipolarons results in a phase transition to a 
superconducting state at a critical temperature in excess 
of a hundred K. Finally, the spin susceptibility and the 
specific heat of the model revealed a separation of charge 
and spin gaps. Because of the presence of a continuum 
spectrum, there is no true ground state gap at any finite 
value of the polaron hopping t. However, strong evidence 
of a finite pseudogap has been found in the range where 
t < J p above T c . 
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